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1. Introduction and Results 



The radiative decay ry' 77 is of special interest in particle theory since it involves 
the gluonic as well as electromagnetic contribution to the axial anomaly. In this paper, we 
show how this decay may be used to measure the topological susceptibility x(0) in QCD 
with dynamical quarks and thus open a window on the topology of the gluon field. 

The flavour non-singlet decay tt^ — 77 is of course well understood and has played an 
important role in the development of the standard model, providing early evidence for the 
existence of colour. Since the pion is a (pseudo) Goldstone boson for the spontaneously 
broken chiral symmetry of QCD, the decay is easily calculated from the electromagnetic 
contribution to the anomaly in the corresponding axial current. 

However, the theory underlying the flavour singlet decay is less developed and indeed 
the special new features arising from the gluon contribution to the divergence of the ^7^(1) 
axial current are usually still ignored in phenomenological analyses. In a previous paper fll], 
we presented an analysis of the theory of rj' —>■ 77 decay in the chiral limit of QCD, 
taking into account the gluonic anomaly and the associated anomalous scaling implied by 
the renormalisation group. Here, we extend that analysis to QCD with massive quarks, 
incorporating 77 — 77' mixing. In particular, we show how a combination of the radiative 
decay formula and a generalisation of the Witten-Veneziano mass formula[|],^ for the 77' 
can be used, under reasonable assumptions, to measure the gluon topological susceptibility 
x(0) in full QCD with massive quarks. 

Our main result is summarised in the following two formulae: 

ex 

9v"i-y + 2nfAgGj-ySao = a^^- (1.1) 

TV 

which describes the radiative decays, and 

(mu{uu) ^ \ 

md{dd) +{2nffA5ao5w (1-2) 

ms{ss) ) 

which defines the decay constants appearing in ( |1 . 1|) through a modification of Dashen's 
formula to include the gluon contribution to the Ua{^) anomaly. 

In these formulae, 77" denotes the neutral pseudoscalars , rj, 77'. The (diagonal) mass 
matrix is {m?)ot/3 and fi'T^^^-y is the appropriate coupling, defined as usual from the decay 
amplitude by 

(77IO = fr?"77 ^\pc.pPlP2^^{Pl)^^{P2) (1-3) 

in obvious notation. The constant a'^^ is the coefficient of the electromagnetic contribution 
to the axial current anomaly: 

d'^J^, = dacbm^QToT'q + 2nf5ao^trG''''G^, + at^^F^^F^^ (1.4) 

OTT OTT 
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where a = 0, 3, 8 is the flavour index and the (i-symbols are deflned from the anticom- 
mutation relations of the generators, {T",T^} = dabcT'^- G^iu and F^^ are the gluon and 
photon fleld strengths respectively. (More precise deflnitions and further notation are given 
in sect. 2.) 



The decay constants in (|1.1| ) are deflned by the relation (|1.2| ) . Notice immediately 
that the decay constants which enter the formula ( |1 . 1| ) are not deflned by the coupling of 
the pseudoscalar mesons tt*', ?7, rj' to the axial current [|T|. In the flavour singlet sector, such 
a deflnition would give a RG non-invariant decay constant which would not coincide with 
the quantities arising in the correct decay formula (|1.1| ) . In practice, since flavour SU{2) 
symmetry is almost exact, the relations for tt*^ decouple and are simply the standard ones 
with /"^^ identifled as /^r, viz. 

O TT 

together with the Dashen formula 

/^m^ = -{mu{uu) + md{dd)) (1.6) 

In the octet-singlet sector, however, there is mixing and the decay constants form a 2 x 2 
matrix: 

The four components are independent. In particular, for broken SU{3), there is no reason 
to express as a diagonal matrix times an orthogonal rj — rj' mixing matrix, which 
would give just three parameters. Several convenient parametrisations may be made, 
e.g. involving two constants and two mixing angles|Q,^^, but this does not seem to reflect 
any special dynamics. 



The novelty of our results of course lies in the extra terms arising in ( p. . 1|) and (1^) due 



to the gluonic contribution to the Ua{^) anomaly. The coefficient A is the non-perturbative 
number which specifies the topological susceptibility in full QCD with massive dynamical 
quarks. Defining the topological susceptibility as 



X(0) = j d'^x mT Q{x) Qm^) {l.i 



where Q = ^trG'^'^ G is the gluon topological charge, the anomalous chiral Ward identi- 
ties determine its dependence on the quark masses and condensates up to an undetermined 
parameter, viz.0 

-A mumdms{uu){dd){ss) 

mumdms{uu){dd){ss) — A{mumd{uu){dd) + mums{uu){ss) + mdms{dd) {ss)) 

(1.9) 
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or in more compact notation, 

Notice how this satisfies the weU-known result that x(0) vanishes if any quark mass is set 
to zero. 

The modified Dashen formula is in fact a generalisation of the Witten-Veneziano 
mass formula for the 77'. Here, however, we do not impose the leading order in 1/Nc 
approximation that produces the Witten-Veneziano formula. Recall that, for nj = 3 and 
non-zero quark masses, this states 

mj, +mj -2m^ = -^x(0)|ym (l-H) 

J TT 



To recover ( p. .Ill ) from our result (see the first of eqs(1.13)) the condensate ms(ss) is 
replaced by the term proportional to f^m^ using a standard Dashen equation, and the 
singlet decay constants are set to ^Jlrifj^^. The identification of the large limit of the 
coefficient A with the non-zero topological susceptibility of pure Yang-Mills theory may 
be seen in different ways, either from the large Nc counting rules quoted below or perhaps 
most simply in the effective Lagrangian analysis explained in section 3. 

The final element in ( |1 . IJ ) is the extra 'coupling' go-i-i in the fiavour singlet decay 
formula, which arises because even in the chiral limit the 77' is not a Goldstone boson 
because of the gluonic U anomaly. A priori, this is not a physical coupling, although 
(suitably normalised) it could be modelled as the coupling of the lightest predominantly 
glueball state mixing with 77'. However, this interpretation would probably stretch the basic 
dynamical assumptionsS underlying ( |1 . 1|) too far, and is not necessary either in deriving 
or interpreting the formula. In fact, the gc-i-i term arises simply because in addition to 
the electromagnetic anomaly the divergence of the axial current contains both the quark 
bilinear operators (^5 = (pj^T^-q and the gluonic anomaly Q. Diagonalising the propagator 



^ The (standard) dynamical assumptions made in deriving (|1.1| ) are described more com- 



pletely in the following sections. Essentially, (LI) is based on the zero-momentum anomalous 
chiral Ward identities. It is then assumed that the decay 'constants' f°'°'{k'^) and couplings 
Qri'^-yyik^) arc approximately constant functions of momentum in the range from zero (where the 
Ward identities are applied) to the relevant physical mass. Notice that this is only applied to pole- 
free quantities which depend only implicitly on the quark masses. This is of course simply the 
standard PCAC or chiral Lagrangian assumption, and corresponds to assuming pole-dominance 
of the propagators for the appropriate operators by the pseudo Goldstone bosons. It becomes 
exact in the chiral limit. 
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matrix for these operators isolates the ry and rl poles, whose couplings to 77 give the usual 
terms g^q-^-^ and Qr^'-yy However, the remaining operator (which we call G - see following 
sections) also couples to 77 and therefore also contributes to the decay formula, whether 
or not we assume that its propagator is dominated by a 'glueball' pole. 



Of course, the presence of the in general unmeasurable coupling go-y-y in (|1-1| ) appears 
to remove any predictivity from the rj' —>■ 77 decay formula. In a strict sense this is true, 
but we shall argue below that it may nevertheless be a good dynamical approximation 



to assume go-y-y is small compared to (7,7'7-y. In this case, we can combine eqs.(|LT|) and 
( p..2|) to give a measurement of the non-perturbative coefficient A in x(0). Assume flavour 
SU{2) is exact so that eqs. (|1.6| ) and ( |1.7| ) for the pion decouple. The remaining parts of 



( |1 . 1|) and (1^) then together provide five equations (since (1^) is symmetric), in which we 



assume the physical quantities m^, m^', grj-y-y and gri' -y-y are all known and we neglect ga-yy- 
Four of these equations may be used to determine the four decay constants The final 
equation is the flavour singlet Dashen formula, which may then be solved for A. This is 
the generalisation of the Witten-Veneziano formula. 

Without neglecting go-y-y, the flve equations give a self-consistent description of the 
radiative decays, but are non-predictive. It is therefore important to analyse more carefully 
whether it is really legitimate to neglect (/c-y-y.il The argument is based on the fact that 
go-y-y is both OZI suppressed and renormalisation group (RG) invariant[Q. Since violations 
of the OZI ru iJ are associated with the Ua{^) anomaly, it is a plausible conjecture that 
we can identify OZI- violating quantities by their dependence on the anomalous dimension 
associated with the non-trivial renormalisation of J^^ due to the anomaly. In this way, 
RG non-invariance can be used as a flag to indicate those quantities expected to show 
large OZI violations. If this conjecture is correct, then we would expect the OZI rule 
to be reasonably good for the RG invariant go-y-y, which would therefore be suppressed 
relative to (7r7'77-0 (An important exception is of course the i]' mass itself, which although 
obviously RG invariant is not zero in the chiral limit as it would be in the OZI limit 



^ Notice that this approximation is implicitly made in almost all the phenomenological 



analyses of radiative pseudoscalar decays, which use a formula like (1.1) but omitting the gc-yy 
term and (mistakenly) implying that the decay constants are defined as the couplings of the axial 
current to 77 and rj' . 

^ The OZI approximation to QCD may be given a precise definition as the truncation of 
full QCD in which non-planar and quark-loop diagrams are retained, but the diagrams which give 
purely gluonic intermediate states (those in which the external currents are attached to different 
quark loops) are omitted. This is a closer approximation to QCD than, for example, the leading 
large Nc limit. 

^ On the other hand, we would not expect the RG non-invariant and anomaly sensitive 
'decay constant' f^^ defined by {0\J%\ti') = ik^f^^ to be well-approximated by its OZI value. 
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of QCD.) Notice that this conjecture has been apphed aheady with some success to the 
'proton spin' problem in polarised deep inelastic scattering P, p!0| , p!T| . 

It is also interesting to look at eqs. (|l.l| ) and ( |1.2D from the point of view of the large 
Nc expansion. The large Nc counting for the various quantities involved is as follows: 
r° = OiVK). = OiVK). 9Gj-y = m2, = 0(l/iV,), ml = 0(1), {qq) = 

0{Nc), Ogjjj = 0{Nc) and A = 0(1). Notice first that this implies x(0) ^ A in the large 
Nc limit, as already used above to derive the Witten-Veneziano formula. It also follows 
that the term in the decay formula involving Ago^^ is suppressed by 0{1/Nc) relative to 
the /""5'r;'^77 term. If large A''^ is reliable in this case, we would therefore indeed expect 
this contribution to be suppressed. However, the large Nc limit must be used with great 
caution in the channel. For example, the same argument would equally imply that 

the additional term A in the flavour singlet Dashen formula is suppressed by 0{1/Nc), yet 
we know that although formally of 0{1/Nc), m^, is not small phenomenologically. 

To make the phenomenological application of our results (|1.1| ) and (|1.2| ) quite clear, 
we now write out the five equations in the 77 — 77' sector explicitly. Set 71/ = 3 and take 
niu = rud = for simplicity. The decay equations are: 



' ^ 

where a'^^ = |A^c and a^-^ = -^^-^c and the Dashen equations are: 



" TT 

' "em „ 
TT 



(1.12) 



{f'^'Yml, + {f^Yml = -4m,{ss) + 36 A 
r^'r^ mj, + rV'X = (1.13) 

Clearly, the two purely octet formulae can be used to find f^"^ and f^"^ if both grj-y-y and 
grf'-y-y are known. The off-diagonal Dashen formula then expresses f^"^ in terms of f'^^ . This 
leaves the two purely singlet formulae involving the still-undetermined decay constant f^^ , 
the topological susceptibility coefficient A, and the coupling gc-yy- The advertised result 
follows immediately. If we neglect go-y-y, we can find /^^ from the singlet decay formula 
and thus determine A from the remaining, generalised Witten-Veneziano, formula. 

On the other hand, we may regard A as a number to be predicted by non-perturbative 
theoretical calculations. In that case, the Dashen formula determines the decay constant 
f^"^ in terms of A, in which case everything is known in the singlet decay formula except 
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the coupling (/c-y-y, which is therefore predicted. Unfortunately, it is not at all clear how 
this could be compared to an experimental measurement without invoking 'glueball dom- 
inance' in the {G G) propagator and assuming that this dominates over other neglected 
pseudoscalar poles in the propagator matrix (see section 3). 

Finally, we should discuss briefly whether it is possible to determine any of the quan- 
tities or A by non-perturbative calculations, either using QCD spectral sum rules or 
lattice gauge theory. The decay constant definition ( p..2| ) has an alternative (prior) form 
in terms of the propagators for the pseudoscalar quark bilinears cj)'i^. In fact (for notation 
and the derivation see sections 2 and 3), we have: 

n{m^)apf^^' = dacem (</>5 (p^)^ ddBf{<pf) (1.14) 

where on the r.h.s. (^5 denotes the (cd) component of the inverse of the matrix 

of two-point functions of the pseudoscalars, taken at zero momentum, and are the 
usual quark condensates. Comparing with ( |1.2[ ) shows that a successful calculation of the 



r.h.s. would imply a determination of the coefficient A governing the topological suscepti- 
bility. Although this looks relatively straightforward, perhaps unsurprisingly it turns out 
to be a very delicate calculation indeed in the QCD spectral sum rule approach,^ primarily 
because the effects of gluons and the anomaly have to make an important contribution to 
what are in first approximation purely quark bilinear propagators. This may nevertheless 
be an interesting problem to pursue or to study on the lattice. 

The other possibility is to determine A by a direct calculation of the topological 
susceptibility x(0) in full QCD with dynamical, massive quarksil. Because of the intricate 
dependence on all the quark masses (including the light quarks) this does not seem feasible 
using spectral sum rules. The situation may be better on the lattice, however, if the 
possibility to change quark masses is exploited. Perhaps A could be extracted from a 
calculation of x(0) in the limit of equal, but not too light, quark masses. 

A general point highlighted once again by this is that it would be extremely useful to 
discover a calculational method that yields IPI vertices directly rather than deducing them 
by amputation of the related Green functions. For example, as we show in section 2, the 
coefficent A in the topological susceptibility is in fact just given by the 'two-point vertex' 
functional Tqq. These 'IPI' functionals are smooth and free of all the singularities and 



^ I would like to thank Stephan Narison for collaborating on a preliminary investigation of 
this problem. 

" bee refs.|l|,|3| for reviews of progress in calculating the topological susceptibility in full 



QCD on the lattice. Ref . [p^ contains results and references relevant to the spectral sum rule 
approach. (See also ||To|JTl|] .) 
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delicate mass dependence associated with external propogator poles. They are therefore 
the essential non-perturbative quantities we need to find. 



It would also be interesting to make a detailed comparison of the formulae presented 
here with the corresponding results in the recently constructed chiral Lagrangians in which 
the 77' is incorporated in the framework of the 1/Nc expansion P,p!5| . Although the for- 



malisms look rather different, mainly because of the non-linear realisations used in the 
chiral Lagrangian approach and the explicit reliance on the l/N^ expansion, the essential 
dynamical input and assumptions are the same and the final physical predictions should 
agree within the limitations of the approximations. 



The rest of this paper is devoted to the derivation of the above results by different 
methods. In section 2 we establish some convenient notation and review the basic anoma- 
lous chiral Ward identities which are the basis of all the subsequent work. Section 3 gives 
a derivation of (|1.1|) and ( |L^ ) using the method of IPI vertex functionals used to analyse 
rj' —>■ 77 in the chiral limit in refs|]l|. Then in section 4 we write an effective Lagrangian 
(extending refs.|^,^) which incorporates all the constraints of the zero-momentum chi- 
ral Ward identities and includes the coupling to electromagnetism. In section 5 we give 
a slightly simplified derivation which follows as closely as possible the traditional PCAC 
methods, generalised as necessary to take account of the gluonic Ua{^) anomaly. This 
section is intended to be as self-contained as possible, and readers interested simply in the 
phenomenology may prefer to go directly to section 5. 



Finally, the methods of this paper may equally be applied to other decays involving 
the r]' [17,^,18|. Of special interest are, for example, r]' V'y where y is a light 1~ meson 
such as p, which is related to 77' —>■ 77 by vector meson dominance; rj' — > 7r7r7, which is 
determined by the box anomaly for one vector and three axial currents; and i/j —>■ r]''j |1T9(] . 
Some interesting current experimental studies of r]' physics include, for example, B rj' K 
by the CLEO collaboration[|g; ip2S il by the BES collaboration at BEPC[|l|; and ry' 
photoproduction at ELSAP[ and CEBAFI 



2. Chiral Ward Identities 



The anomalous chiral Ward identities for QCD with massive quarks have been written 
down in the form used here in ref.|jTlI] and reviewed in [^^. We refer to these papers for 
more complete derivations and in this section simply define our notation and quote the 
essential identities. In this section, we omit the electromagnetic contributions. 

The composite operators involved in the Green functions and IPI vertices studied here 
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are the currents and pseudoscalar operators J^g, Q, (^5 and the scalar (f)"- where 



a. 










rrid 











J^B = qi^l^Ty Qb = ^trG^.G^'^ 

G^i, is the gluon field strength. In this notation, = are fiavour SU{nf) generators, 
and we include the singlet generator = 1 and let the index a = 0, i. We only 

need to consider fields where i corresponds to a generator in the Cartan sub-algebra, so 
that a = 0, 3, 8 for 77,/ = 3 quark flavours, (i-symbols are defined by {T",T''} = dabcT'^- 
Since this includes the fiavour singlet Ua{^) generator, they are only symmetric on the 
first two indices. For = 3, the explicit values are (iooo = <^033 = <^088 = 2, dsso = dsso = 

1/3, dsss = ^383 = —'^888 = 1/ V^- 

We also use the following compact notation. The quark mass matrix is written as 
m"'T°', so that for n/ = 3, 



m^l + m^T^ + m^T^ (2.2) 



In the same way, the chiral symmetry breaking condensates may be written as 

(uu) \ 1 
{dd) = -((/)0)1 + 2((^3^T3 + 2((/)S)T8 (2.3) 
{ss) J ^ 

where (cp'^) is the VEV (qT'^q). We then define 

Mab = dachm'' ^ab = dabc{4>^) (2.4) 

Eq ( |2.1D defines the bare operators. The renormalised composite operators are defined 
as follows: 

tO _ 7 tO ja^O _ ra/O 

Q = QB-^^{l-Z)d^J% (2.5) 

<t)i = z^<t^iB r = z^^B 

where Z^ is the inverse of the mass renormalisation, = Z^^. The non-trivial renormal- 
isation of J^^ means that its matrix elements scale with an anomalous dimension 7 related 
to Z. This occurs because J^^ is not a conserved current, due to the anomaly Q. Notice 
also the mixing of the operator Q with d^J^^^ under renormalisation. 
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The Green functions for these operators are constructed by functional differentiation 
from the generating functional W[V^^,9, S^, S""], where V^^, 9, S^, S"- are the sources 
for the composite operators J^5, Q, (/)5, respectively. For example, the Green function 
i{0\T Q{x) Q{y)\0) is given by se{x)Y8{y) ' '^hich we abbreviate as W00. In this notation, 
the anomalous zero-momentum chiral Ward identities are: 



2nfSaoWee + MacWs^e = 
2nf5aoWQsb + MacWs^s^ + = 



(2.6) 



which implies the following identity for the topological susceptibility, 

(2n/)2x(0) = MocVTgcg.Mrfo + (M$)oo (2.7) 
These are derived from the fundamental anomalous Ward identity 

d^Wv^^ - 2nf5aoWe - M^cWg^^ + dadcS'^Ws^ - dadcStWs^ = (2.8) 

which is the precise expression in the functional formalism of the identity ( |1.4| ) . 

The IPI vertices used in section 3 are defined as functional derivatives of a second 
generating functional (effective action) F, constructed from by a partial Legendre trans- 
form with respect to the fields Q, (^5 and only [not the currents ^^5)@- The resulting 
vertices are 'IFF w.r.t. the propagators for these composite operators only. This separates 
off the particle poles in these propagators, and gives the closest identification of the field 
theoretic vertices with the physical couplings such as Qn'^-yj. 

The basic anomalous chiral Ward identity for F follows immediately from ( p.8|) for W: 

d^Fv^a^ - 2nfdaoQ - Mac(t>l + dacd4^''T^. - dacd(t)iT^. = (2.9) 

and other identities follow simply by functional differentiation. In particular, for the two- 
point vertices, we find the following zero-momentum identities (analogous to 



^acr^cQ - 2nf5aO =0 



which together imply 

'^acT^i^i'^db = -{M^)ab (2.11) 



These will be useful in section 4. 



The fact that the topological susceptibility is zero for vanishing quark mass can be 
seen immediately from (|2.7| ) . One of the simplest ways to derive the precise form ( p..9| ) or 
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( 1.10|) is in fact to use an identity involving F. As is well-known, the two-point vertices are 
simply the inverse of the two-point Green functions (propagators), so in the pseudoscalar 
sector we have the following matrix inversion formula: 



(2.12) 



= -(Wee - Wes^Mac{MWs,s,M)~^MdbWsie] 
and using the identities (|2.6|) and (|2.7| ) this implies 

Tqq = - (2n^)2x(M$)oV)"' (2.13) 

all at zero momentum. Inverting this relation gives the result for the topological suscepti- 
bility: 

X = -rQ^(l - (2n/)2rQ^(M$)o-oi)"' (2.14) 

Substituting the explicit expression for (M$)qq^ (which is easily found from the definitions 
above), viz. 

we see that ( |2.14D reproduces the general form (|1.10| ) where we can now identify the (mass- 
independent) non-perturbative coefficient as 

^ = rQ^ (2.16) 



3. 1]' — > 77 from IPI Vertices 

In this section, we present the most theoretically complete derivation of the decay 
formula ( |1 . 1|) and generalised Dashen formula (|1.2| ) . This follows the derivation previously 
given in ref.[|l| for the chiral limit, extended to include quark masses. The technique relies 
on the identification of the couplings Qn'^-y-y with the zero-momentum limit of certain IPI 
vertex functions, precisely defined using the Legendre transform F introduced in section 
2. These techniques have also been used in our series of papers on the t^A(l) Goldberger- 
Treiman relation and the 'proton spin' P JTO| , pT| . See also p5|j2^ for reviews. 



The starting point is the Ward identity ( p.9| ) extended to include the electromagnetic 
contribution to the anomaly for the axial current: 

d^Tv^a^ - 2nf5aoQ - <^Qem(A) - Mac(t)l + f^acd'^'^F^c - cf„,rf(^^r^c = (3.1) 
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where QemiA) is shorthand notation for -^Ff^^F'^'^ , where is the field strength for 
the electromagnetic field A^. (Since we are working only to leading order in a, it is 
not necessary to consider Qem as an independent composite operator with non-trivial 
renormalisation. ) 

DiflFerentiating twice w.r.t. the field A^, evaluating at the VEVs, and taking the Fourier 
transform, we find 

a ^ a c 

ikuTya A\y^p + a'^T^—exfjaBPiPo + dabc4>^^(hbA^Ap = (3-2) 

where pi,p2 are the momenta of the photons. To simplify notation, it will be conve- 
nient from now on to define vertices F with the kinematical factors removed, in particular 
T^a^xj^p = —T,j,<^A>'Ap^Xa-apPiP2- Noticc that the mass term in (|3.1D does not contribute 
explicitly to this formula. From its definition as IPI w.r.t. the pseudoscalar fields, the 
vertex Tya^A^^p has no pole at /c^ = (even in the chiral limit) so the first term vanishes 
at zero momentum k, leaving simply 



fc=0 TT 



The first step in converting ( |3.3|) to the decay formula ( |1.1| ) is to identify the physical 
states 77". These appear as poles in the propagator matrix for the four pseudoscalar 
operators Q, ^5 (a = 0, 3, 8). To isolate these poles, we diagonalise the propagator matrix 
in this sector then normalise the three operators coupling to the physical states. 

We therefore define the operator 



(3.4) 



so that by construction the propagators {G (^5) all vanish. (Notice that integrations over 
repeated spacetime arguments are implied in this condensed notation.) Then define oper- 
ators 



rf = C 



(3.5) 



such that the propagator matrix 





V 



\ 





(3.6) 



where 5"° are the sources for the operators rf 
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This change of variable affects the partial functional derivatives in r^a^A^^p in 
which involves -^rs: at fixed Q. In terms of the new variables G, ry" we have 



5r7" 5 6G 6 
+ 



(3.7) 



The decay formula therefore becomes 



The decay constants are identified as 

n = $a6C^'" (3.9) 

In terms of the propagators, we can write (from (|3.6| ) ) 

niWs.sXlf^' = ^aciWs.sJ-J^db (3.10) 
and so at zero momentum 

nmlpf^f"' = ^aciWs,s,U^db (3.11) 

as quoted in ( |1.14|) . 

The remaining steps in finding ( p..l|) and ( |1.2|) are an exercise in manipulating the 



zero-momentum Ward identities ( |2.6| ) . First note that combining the two identities in 
(ill) gives 

MacWs^s^Mdb = -{M<5)ab + {2nf)\{0)SaoSbo (3.12) 



whose a, 6 = component is just (|2.7| ) . Note that {M^)ab is symmetric. Also define 
loo = '^ao'^bo. Then we can write 



^ab{Ws.^s.XbWs^,0 = {W)a4MWs.^s,M)jM,,Wsie 

= -2n/(M$)„J-(M$) + (2n/)2x(0)loo)"'x(0) 

V / cO 

= 2n/x(0)(l - (2nj)2x(0)(M$)oo')"'5„o 
12 



(3.13) 



where in the final step we have used the identification ( p.l3|) . Similarly, 

^aciWs.s.V^dt = i^M),,{MWs.^s,M)-J{M^)at 

= (M$)„J-(M$) + (2n/)2x(0)loo)"'(M$)rf, (3.14) 

V / cd 

= -(M$)„, + (2n/)2rQ^ SaoSw 

This establishes the required results. Substituting ( p.l3|) and ( p.l4|) into ( p?8|) and 
( p.ll| ) we find the decay formula 

^abC^^"" t^^^xj^p + 2nfTQQ Tqj^xj^p = o-tm— (3.15) 

where the decay constants satisfy 

nml^f^' = -{M^)aB + (2n;)2rQ^ SaoSw (3.16) 
The final step is to identify the IPI vertices with the couplings defined in section 1, 

viz. 

It is at this point that the central dynamical assumption is made. In fact, eqs ( p.l5|) and 
( pi.l6| ) are exact identities, following simply from the definitions and the zero-momentum 
chiral Ward identities. To make contact with the radiative decays of the physical particles, 
we must assume in particular that the IPI vertex evaluated at k = accurately approxi- 
mates the physical coupling, which is defined on mass-shelM This requires that t^aj^xj^p 
has only a weak momentum dependence in the range < < m^a. This is reasonable, 
since it is defined to be a pole-free, amputated dynamical quantity. However, as in stan- 
dard PCAC, the assumption is expected to be excellent for the tt but progressively worse 
as the mass of the pseudo-Goldstone bosons increases. The hope here, in common with 
all attempts to include the r]' in the framework of PCAC (including chiral Lagrangians 
with 1/A^c effects included B,|T5|, is that the approximation remains sufficiently good at the 
mass of the 77'. 

It is also important to determine the behaviour of all the quantities appearing in 
these formulae under the renormalisation group. Recall that RG behaviour was a key 
factor in the conjecture that go-y-y may be neglected in first approximation in the decay 



The assumption that the IPI vertices as defined here can be identified with the decay 
coupUngs of the physical particles at all rests on the assumption that the dominant particle poles 



in the pseudoscalar propagator matrix are indeed those of the rj" (see ( ^.61) ) 
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formula (|1 . 1| ) . All the required formulae are given in ref. . The result is that all the 
quantities appearing in the final formulae , (|13) (or alternatively ( |3.15| ) , ( p.l6| ) ) 

are RG invariant. However, notice that this is only true for the IPI vertices evaluated at 
k = (and in fact also on-shell), not for arbitray momenta. The proof is quite intricate, 
but since everything can be read off from ref.[^ (see also |T^) we will not give any further 
explanation here. 



4. Effective Action 



The results in the previous section can be summarised by writing an explicit form 
for the effective action r[Q,(/)5] compatible with all the zero-momentum anomalous chiral 
Ward identities. Of course this adds no new physics, but allows the essential results to 
be read off in a perhaps simpler and more systematic way. The resulting effective action 
is essentially identical to the di Vecchia-Veneziano0, Rosenzweig-Schechter-Trahern |1T6[] 
Lagrangian (with ^ = 0) though without explicit reference to the 1/Nc expansion. Less 
obviously, it is also very closely related to the (non-linear) chiral Lagrangians which incor- 
porate the rj' in the framework of large Nc P,P^ . 

I) , (S) has 



The simplest effective action r[Q, (p"^] compatible with the identities 

It is 



been written down in ref. [11 



r[Q,€] 



dx 



1 

2A 



(4.1) 



The final term is written in matrix notation. and fx^^p are matrices, where fx^^p 

defined by the D ashen formula 



IS 



cTi3b 



(4.2) 



The decay constants will subsequently be identified with those in section 3 (as will the 
constant A), is of course the pseudo-Goldstone boson mass matrix in the OZI limit of 
QCD, i.e. before including the coupling to the gluonic anomaly operator Q. 



In (|4.1|) the simplest choice of kinetic terms for the fields (^5 has been made, with the 
chosen to be constants. This is where the dynamical, pole-dominance, assumption of 
standard PCAC (or chiral Lagrangians) is built in. No kinetic terms are included for the 
composite operator Q (no glueball poles), and no higher order terms in Q are included 
(these would be suppressed for large Nc)- 



14 



We have also included terms in ([4.1j ) involving Qem{A) to satisfy the anomalous chiral 
Ward identities with the additional electromagnetic contribution. The term coupling Qem 
to is required, whereas the term QQeui{A) is permitted. (Like A, B is taken to be a 
constant.) Differentiating ( [4.1|) , we immediately obtain 

a. 

^ab r^^b^A^p = al^- (4.3) 



as previously found in 

The second derivatives oiV[Q, (j)"}] are 



(4.4) 



which clearly satisfy the identities ( |3.9| ) and ( p.ll| ) . We confirm the identification Tqq = 
A. The corresponding Green functions are found by inversion: 

Wee = -AA-i 
Wes^^ = 2nfA A'^ 

/ (4-5) 

Ws^e = 2nfA^ac{f{k''-^i^)f] A'^ 

V / cO 

Ws^Sl = -$ac A J ^db 



where 



and 



A = l-{2nf)'A[f{e-i,')f)~'^ (4.6) 



A = f{k'-fi')f-i2nf)^Aloo (4.7) 

However, in this form the propagator matrix is clearly not diagonal and the operators 
are not normalised so as to couple with unit decay constants to the physical states. It is 
therefore convenient to make a change of variables in F so that it is written in terms of 
operators which are more closely identified with the physical states. Of course this change 
of variables is precisely that described already in section 3. We therefore define 

G = Q-Wes^,{Ws,sXb<Pl 
= Q + 2nfA%,'(Pl 

using (|4.5|) for the propagators (c.f. the identity (|3.13D ), and 



^ fT^a^-l^l (4.9) 
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In terms of these operators, the effective action is 
1 



dx 



2A 



1 

H — V 
2 ' 



{-d'-lx')-{2nffAf-H,,r' 



(4.10) 



It is then straightforward to read off the propagators 

{G G) = ~A 



-1 



5"^ 



(4.11) 



where the diagonal mass matrix for the physical 77" states satisfies the generalised Dashen 
fomula 

nmlpfP' = niilpf^' + {2nffA5ao5m (4.12) 

Of course this is identical to ( |1 . 2\j . This confirms the identification of the decay constants 
in the effective action with those in sections 1 and 3. 

Finally, to obtain the decay formula itself, we take functional derivatives of F to get 



OL 

/"'"^■n'^A^AP = Km 2nfABdaQ 

TV 

Combining these, we find 

r^a^A^p + 2nfA Fq^a^p 



(4.13) 



a ^ 

n — 

em _ 
TT 



(4.14) 



in agreement with ( |1 . 1| ) . Notice the importance of including the QQem coupling in 
in obtaining this result. 



5. Ua{1) PCAC 

For the third of our variations on a theme, we present a derivation of the decay 
formula following as closely as possible the traditional language of PCAC. This should be 
reasonably self-contained, though we will use the compact notation defined at the start of 
section 2 and simply quote the chiral Ward identities without proof. 

Consider first QCD by itself without the coupling to electromagnetism. The axial 
anomaly equation is 

d''J^,=Mab(pt + 2nfQ5ao (5.1) 
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where J^g is the axial current, (^5 the pseudoscalar quark bihnear operator and Q the 
topological charge. Mab describes the quark masses and the condensates. The anoma- 
lous chiral Ward identities, at zero momentum, for the propagators (i.e. two-point Green 
functions) of these operators are 



2n;(Q Q)daQ + Macm Q) = 
2nf{Q (Pl)6ao + Mac{(pt (Pl) + ^ab = 



(5.2) 



which imply 

MacMbdm (pi) = -{M^)at + {2nff{Q Q)5aQ5to (5.3) 
We also need the result for the general form of the topological susceptibility: 

X(0) = {Q Q) = — ^ (5.4) 

^ 1 - (2n/)2A(M$)ooi ^ ^ 

Although the pseudoscalar operators (^5 and Q indeed couple to the physical states 
77" = ry', 77, tt'^, it is more convenient to redefine linear combinations such that the resulting 
propagator matrix is diagonal and properly normalised. That is, we define operators ry" 
and G such that 

{QQ) {Q^l)\ ^ ({GG) 



(5.5) 



This is achieved by 



and 



G = Q-{Q ct>t){{ct>, ct>,))-,'4 



(5.6) 



^ fT^a^-l^l (5.7) 

With this choice, the {G G) propagator is 

{G G) = -A (5.8) 

and we impose the normalisation 



This implies that the constants in ( ^.7|) , which we see shortly are simply the decay 
constants, must satisfy the (Dashen) identity 

nmlpff"' = ^aM, <P,)V^db 

= -{M^)ab + i2nf)^A5aoSbo 
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The last line follows from the Ward identities (|5.3|) and ( |5.4|) . In terms of these new 



operators, the anomaly equation (|5.1| ) now reads simply: 

d^J^ = nKpv" + 2nfGdao (5.11) 

After these preliminaries, we now recall how conventional PCAC is applied to the 
calculation of tt^ 77. The pion decay constant is defined as the coupling of the pion to 
the axial current 

{0\j'^,\r:) = ikj^ {0\d''jf,,\7r) = Uml (5.12) 

and satisfies the Dashen formula 

f^i^l = -{rnu + md){qq) (5.13) 

The next step is to define a 'phenomenological pion field' tt by 

d^J% ^ UmlTX (5.14) 

This is the step at which the crucial 'pole-dominance' assumption is made. Now include 
electromagnetism. The full anomaly equation is extended as in ( |1.4| ) to include the F^'^F^i, 
contribution. Using ( |5.14| ) we therefore have 

a ~ (5-15) 

= Umlin 7r)(77K) + <^-(77|i^^'^i^M.|0) 

where (tt tt) is the pion propagator — 1/(A;^ — m^). At zero momentum, the l.h.s. vanishes 
because of the explicit factor and the absence of massless poles. We therefore find, 
defining the couplings as in ( p. .31 ) , 

TT 

In the full theory including the fiavour singlet sector and the gluonic anomaly, we 
find a similar result. The 'phenomenological fields' are defined by (|5.11| ) where the decay 



constants satisfy the generalised Dashen formula (|5.10|) . Notice, however, that they are not 



simply related to the couplings to the axial current as in ( |5.12|) for the fiavour non-singlet. 
We therefore find: 



(5.17) 
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using the fact that the propagators are diagonal in the basis ry" , G. Using the exphcit 
expressions ( p.8| ) and ( pT9| ) for the propagators, evaluating at zero momentum, and setting 
the l.h.s. to zero, we find in this case: 

a 

f"" g^a^^ + 2nfA gc-r-y SaO = «em- (5.18) 

where the extra coupling ga-y-y is defined through ( ^.ITp . This completes the '^7^(1) 
PCAC derivation. It is evidently a straightforward generalisation of conventional PCAC 
with the necessary modification of the usual formulae to take account of the extra gluonic 
contribution to the axial anomaly in the fiavour singlet channel, the key point being the 
identification of the operators ry" and G in ( ^.11|) . 
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